Generic Miniband Structure of Graphene on a Hexagonal Substrate 
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We present a symmetry based classification of generic miniband structures for electrons in 
graphene placed on substrates with the hexagonal Bravais symmetry. In particular, we identify 
conditions at which the first moire miniband is separated from the rest of the spectrum by either 
one or a group of three isolated mini Dirac points and is not obscured by dispersion surfaces coming 
from other minibands. In such cases the Hall coefficient exhibits two distinct alternations of its sign 
as a function of charge carrier density. 

PACS numbers: 73.22.Pr,73.21.Cd,73.43.-f 



Recently, it has been demonstrated that the electronic 
quality of graphene-based devices can be dramatically 
improved by placing graphene on an atomically flat crys- 
tal surface such as hexagonal boron nitride (hBN) 
[7]. At the same time, graphene's electronic spectrum 
also becomes modified, acquiring a complex, energy- 
dependent form caused by incommensurability between 
the graphene and substrate crystal lattices |8lu0|. In 
particular, for graphene placed on hBN, the difference 
between their lattice constants and crystallographic mis- 
alignment generate a hexagonal periodic structure known 
as a moire pattern [2 02 I8HTU] . The resulting periodic 
perturbation, usually referred to as a superlattice, acts 
on graphene's charge carriers and leads to multiple mini- 
bands and the generation of secondary Dirac-like spectra. 
The resulting new Dirac fermions present yet another 
case where graphene allows mimicking of QED phenom- 
ena under conditions that cannot be achieved in particle 
physics experiments. In contrast to relativistic particles 
in free space, the properties of secondary Dirac fermions 
in graphene can be affected by a periodic sublattice sym- 
metry breaking and modulation of carbon-carbon hop- 
ping amplitudes, in addition to a simple potential mod- 
ulation. The combination of different modulation mech- 
anisms results in a multiplicity of possible outcomes for 
the moire miniband spectrum in graphene which we sys- 
tematically investigate in this Letter. 

To describe the effect of a substrate on electrons in 
graphene at a distance, d, much larger than the spacing, 
a, between carbon atoms in graphene's honeycomb lat- 
tice, we use the earlier observation IMT41 that, at d 3> a 



the lateral variation of the wavefunctions of the p z car- 
bon orbitals is smooth on the scale of a. This is man- 
ifested in the comparable sizes of the skew and vertical 
hopping in graphite and permits an elegant continuum- 
model description (TTHT3] of the interlayer coupling in 
twisted bilayers and the resulting band structure. A sim- 
ilar idea applied to graphene on a h-BN substrate [8HTU] 
suggests that a substrate perturbation for Dirac electrons 
in graphene can be described in terms of simple harmonic 
functions corresponding to the six smallest reciprocal lat- 
tice vectors of the moire superlattice. 

Below, we shall use a similar approach to analyse the 
generic properties of moire minibands for electrons in 
graphene subjected to a substrate with a hexagonal Bra- 
vais lattice with a slightly different lattice constant of 
(l+6)y/3a, \S\ <C 1, compared to that of y/3a for graphene, 
and a small misalignment angle, 6<1. The moire pat- 
tern harmonics are described by vectors b m = R^m/sb® 
(to = 0,1... 5), which can be obtained by the anticlock- 



wise rotation, 
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with |6o| = fr~ §f VS 2 +9 2 . For a substrate with a simple 
hexagonal lattice or a honeycomb lattice with two iden- 
tical atoms, the perturbation created for graphene elec- 
trons is inversion-symmetric. For a honeycomb substrate 
where one of the atoms would affect graphene electrons 
stronger than the other (e.g. such as h-BN, for which the 
occupancy and size of the p z orbitals are different) the 
moire potential can be modelled as a combination of a 
dominant inversion-symmetric part with the addition of 
a small inversion-asymmetric perturbation, 
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I 



The Hamiltonian, H, acts on four-component wave- functions, (^ak,^bk,^bk',—^ak') T , describing the 
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FIG. 1: (a) The hexagonal Brillouin zone for the moire superlattice. (b) Three volumes in the space of the moire superlattice 
parameters where the edge of the first miniband, in graphene's valence band, contains an isolated mDP at the K-point (red) or 
the — K-point (blue) or three isolated mDPs at the sBZ edge (green). Parameters for which the ±ft-point is triple degenerate 
are shown by the red and blue surfaces. (c)The same for the conduction band in graphene. 



electron amplitudes on graphene sublattices A and B 
and in two principal valleys, K and K' . It is written 
in terms of direct products <JiTj, of Pauli matrices <7j 
and Tj separately acting on sublattice and valley indices. 
The first term in H is the Dirac part, with p= i\7 — eA 
describing the momentum relative to the centre of the 
corresponding valley, with Vxi = B. The rest of the 
first line in Eq. ([!]) describes the inversion-symmetric 
part of the moire perturbation, whereas the second line 
takes into account its inversion-asymmetric part. In the 
first line, the first term, with fi(r) = J2 m =o 5 elbm T i 
describes a simple potential modulation. The second 
term, with /2(r) = iJ2 m =o s( — i) m e ibm " n , accounts for 
the A-B sublattice asymmetry, locally imposed by the 
substrate. Then the third term, with unit vector l z , de- 
scribes the influence of the substrate on the A-B hopping: 
consequently [T5lll7| . this term can be associated with a 
pseudo-magnetic field, e(3 — ±mb 2 f2(r), which has op- 
posite signs in valleys K and K' . All parameters <C 1 
111 Eq. are dimensionless and vb rj 2tt^/5 2 + 2 jo, 
where 70 ~ 3eV is the nearest neighbour hopping inte- 
gral in the Slonczewski- Weiss tight binding model |18j . 
Concerning the inversion-asymmetric part, the second 
line in Eq. ([I]), we assume that \ui\ <§; |itj|. Note that 
the last term in each line can be gauged away using 

In the absence of a magnetic field, the Hamiltonian 
Eq. ([T]) obeys time-reversal symmetry, which follows from 
both (Ji and Tj changing sign upon the transformation 
t — >• — t [12]. As a result, £k+p — ^K'-p and we limit the 
discussion of minibands to the K valley. Moreover, using 
the commutation properties of <ii one can establish that 
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To calculate the miniband spectrum for H in Eq. we 
perform zone folding (in the graphene K valley) bring- 
ing states with momenta related by the reciprocal lat- 



tice vectors nib\ + 712^2 of the moire pattern to the 
same point of the superlattice Brillouin zone (sBZ) in 
Fig.flja). Then, we calculate the matrix elements of H 
between those states and diagonalise the corresponding 
Heisenberg matrix numerically exploring the parametric 
space (uq, ui, U3) of the moire superlattice shown in Fig 
[l](b,c). The size of the matrix is chosen to guarantee 
the convergence of the calculated energies for the three 
lowest minibands in both the conduction band (s = +1) 
and the valence band (s = —1). Below, we discuss the 
generic features of the moire miniband spectra for the 
characteristic points in the parametric space (ito, Mi, 1*3), 
identified by black dots in Fig. [l|b,c), using both the 
numerically calculated dispersion surfaces in Fig. [2] and 
analytical perturbation theory analysis. 



For the zero-energy Dirac point in graphene, there are 
only the original p = states in each valley that appear 
at e = upon zone folding. For all three characteris- 
tic spectra shown in Fig. [2] for the inversion symmetric 
moire perturbation, the gapless Dirac spectrum persists 
at low energies near the conduction-valence band edge 
with almost unchanged Dirac velocity, [l + 0{u 2 )\ v, 
whereas the the inversion asymmetric terms are able [10] 
to open a minigap. 



For the point fi = bo/2 on the edge of the first sBZ, 
zone folding brings together two degenerate plane wave 
states, \fi + q) and \fi + b 3 + q). The splitting of these 
degenerate states by the moire potential in Eq. ([I]) can 
be studied using degenerate perturbation theory. The 
corresponding 2x2 matrix, expanded in small deviation 
q of the electron momentum from each of the three sBZ 
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FIG. 2: Numerically calculated moire miniband (left), the corresponding density of states (centre), and Landau level spectrum 
(right) for electrons in the vicinity of graphene's K point 



/x-points [30] has the form 
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Hu « (sui - u 3 ) - i(sui - u 3 ) + 2y (uq + iu ). (3) 

For the inversion-symmetric perturbation, the dispersion 
relation resulting from Eq. (13) contains an anisotropic 
mini Dirac point (mDP) [5J Oil H2 with Dirac velocity 
component w 2uqv in the direction of the sBZ edge and 



w v in the perpendicular direction. This feature is clearly 
seen at the /i-point of the first moire miniband in the 
valence band, in the top row of Fig. [2] Note that the 
electron spectrum is not symmetric between the valence 
and conduction bands and that the mDPs at the /x-point 
in the conduction band are obscured by an overlapping 
spectral branch. 

Moving in parameter space, e.g., along the line shown 
in Fig. [ijb) , the positions of the three anisotropic mDPs 
shift along the sBZ edge towards the sBZ corners: ei- 
ther k — (64 + b 5 )/3, or —k, as shown by arrowed 
lines in Fig. [lja). In general, a spectrum with three iso- 
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lated mDPs at the sBZ edge is typical for the green vol- 
ume in the parameter space in Fig. [ljb) for the valence 
band, or Fig. ITTc) for the conduction band. In contrast, 
for (uo,Mi,U3jin the clear part of the parameter space, 
mDPs on the edge of the first sBZ are overshadowed by 
an overlapping spectral branch, as is the case on the con- 
duction band side for all three cases shown in Fig. [2] 

For the points in Fig. [ljb,c) on the red and blue sur- 
faces, the three mDPs reach the K-point, forming a triple 
degenerate band crossing, as in the valence band spec- 
trum shown in the middle row of Fig. [2] which can be 
traced using the perturbation theory analysis of the band 
crossing at k discussed below. 

The third line in Fig. [2] shows the third type of spec- 
trum of moire minibands, characteristic for the red and 
blue volumes of the parameter space in Fig. [T] The 
characteristic feature of such spectra consists in a sin- 
gle isolated mDP, at the ±K-point, in the valence band 
(Fig. |2|b)) or the conduction band (Fig. [l|c)). 

For the k and — K-points, zone folding brings together 
three degenerate plane wave states, \((n + bi + 

q)), and |£(k + b 2 + q)) (where £ = ±), whose splitting 
is determined by 
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For 7^ 0, the inversion-symmetric terms in Eq. Q 
partially lift the £/«-point degeneracy into a singlet with 
energy E K — 2wq and a doublet with energies E K + , 
so that a distinctive mDP [5J characterised by Dirac 
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FIG. 3: The relation between the two densities at which the 
Hall coefficient in graphene reverses sign upon it's doping with 
holes. The results are shown for several realisations of moire 
superlattice in the parameter range corresponding to either 
three isolated mDPs on the sBZ edge (squares) or one isolated 
mDP at the sBZ corner (other symbols). The thresholds for 
isolation are indicated on the x-axis. 



velocity v K = [1 + 0(u)] | is always present at ±k some- 
where in the spectrum [53] . Note that each isolated mDP 
is surrounded by Van Hove singularities in the density of 
states corresponding to saddle points in the lowest en- 
ergy minibands, and, that the presence of the weaker 
inversion-asymmetric perturbation is able to open a mini- 
gap in the mDPs. 

Also, the appearance of mDPs at the edge of the first 
miniband results in a peculiar spectrum of electronic 
Landau levels, as shown on the r.h.s of Fig. [2] Each 
data point in these spectra represents one of the Hofs- 
tadter minibands |24j (with an indistinguishably small 
width) calculated for rational values of magnetic flux, 
|$o P er moire supercell following a method in Ref. |13j . 
Using these spectra one can trace a clearly separated 
"zero-energy" Landau level related to the isolated k- 
point mDP in the valence band in the bottom row of 
Fig. [2j in addition to the true zero-energy Landau level 
at the conduction-valence band edge. The three isolated 
mDPs on the sBZ edge in the valence band (top row 
of Fig. [2]) also results in a "zero-energy" Landau level, 
though not as clearly separated and split by the mag- 
netic breakdown occurring at $ s» 0.1<I>o- 

To summarise, the inversion symmetric moire pertur- 
bation will result in either the first sBZ separated from 
the rest of the spectrum by one or three mDPs, or, for 
weak perturbations, will result in overlapping first and 
higher minibands. The experimental consequences of this 
can be expected in the optical spectroscopy of graphene 
on a hexagonal substrate: the presence of mDPs and Van 
Hove singularities in the density of states should lead to 
a modulation of the FIR and IR absorption spectra of 
monolayer graphene, which otherwise, has the flat ab- 
sorption coefficient of 2.3%. Another experimental con- 
sequence of the moire minibands would consist in a non- 
monotonic variation of the Hall coefficient upon doping 
the graphene flake with electrons or holes. For exam- 
ple, for those miniband spectra in Fig. [2] where there are 
isolated mDPs in the valence band, the Hall coefficient 
would pass through a zero value and change sign at two 
characteristic Hall densities, n\ and At the density n\ 
which corresponds to the valence band filled with holes 
up to the Van Hove singularity the Hall coefficient will 
change sign from positive to negative. At the higher Hall 
density, n 2 , which corresponds to a completely filled first 
miniband, it would repeat the behaviour at the neutrality 
point changing sign from negative to positive. Such be- 
haviour is expected to take place for the entire regions of 
the parametric space painted red, blue or green in Fig. [T] 
The relation between these two carrier densities for vari- 
ous types and strengths of moire perturbations is shown 
in Fig. [3j For the clear part of the parametric space for 
which we find substantial overlap between many moire 
minibands such alternations in the sign of the Hall coef- 
ficient would be obscured by the competing contributions 
from the "electron-like" and "hole-like" branches in the 
spectrum. 
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